Random matrix ensembles from nonextensive entropy 
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The classical Gaussian ensembles of random matrices can be constructed by maximizing 
Boltzmann-Gibbs-Shannon's entropy, ^bcs = — / dH[P(H)] ln[P(H)], with suitable constraints. 
Here we construct and analyze random-matrix ensembles arising from the generalized entropy 
Sq = {l- j dH[P(H)]'') /{q - 1) (thus 5*1 = S'bgs). The resulting ensembles are characterized 
by a parameter q measuring the degree of nonextensivity of the entropic form. Making q —* 1 recov- 
ers the Gaussian ensembles. If g 7^ 1, the joint probability distributions -P(H) cannot be factorized, 
i.e., the matrix elements of H are correlated. In the limit of large matrices two different regimes 
are observed. When q < 1, -P(H) has compact support, and the fluctuations tend asymptotically 
to those of the Gaussian ensembles. Anomalies appear for q > 1: Both -P(H) and the marginal 
distributions P{Hij) show power- law tails. Numerical analyses reveal that the nearest-neighbor 
spacing distribution is also long-tailed (not Wigner-Dyson) and, after proper scaling, very close to 
the result for the 2x2 case - a generalization of Wigner's surmise. We discuss connections of these 
"nonextensive" ensembles with other non-Gaussian ones, like the so-called Levy ensembles and those 
arising from soft-confinement. 



PACS numbers: 
mechanics) 



05.45.Mt (quantum chaos), 24.60.Lz (chaos in nuclear systems), 05.20.-y (classical statistical 



I. INTRODUCTION 



The Gaussian ensembles of Random Matrix Theory 
provide the standard statistical description of spectral 
fluctuations in a multiplicity of quantum systems rang- 
ing from nuclei to disordered mesoscopic conductors and 
classically chaotic systems 0, 0, IS S • 

Gaussian ensembles can be obtained from two postu- 
lates: the invariance of the joint distribution probability 
P(H) with respect to changes of bases and the statisti- 
cal independence of matrix elements. An alternative and 
more appealing way of constructing random matrix en- 
sembles uses a maximum entropy principle 0, 0] . One 
constraint is normalization, 



converges) , 



j dHP(H)trF(H) = 1 



(2) 



(the trace ensures rotational invariance). For instance, 
the Gaussian ensemble of real symmetric matrices is ob- 
tained by the simplest choice 



(3) 



/ 



(iHP(H) = 1 



(1) 



The other one has the purpose of confining the spec- 
trum, but is otherwise arbitrary (as long as the integral 
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It has been proven that, in the limit of large matrices, 
and for a strong enough confining potential V, local fluc- 
tuation properties tend to those of the Gaussian case, 
whatever the shape ofV 0; 13 • 

To escape from Gaussian universality one must con- 
sider soft-conflnement potentials ^ J^] , or breaking ro- 
tational invariance. The latter case typically arises when 
matrix elements Hij are independent (non Gaussian) ran- 
dom variables. For instance, Cizeau and Bouchaud con- 
structed anomalous "Levy ensembles" by drawing Hij 
from a long-tailed distribution pTl |. 

The purpose of this paper is to present a new way 
of constructing non-Gaussian ensembles while preserving 
rotational invariance. The idea is to use a maximum 
entropy approach with the usual constraints but with the 
nonextensive entropy |l2j| : 



l-/rfH[P(H)]^ 

q-1 



(4) 



2 



where g is a free real parameter (q = 1 recovers Shannon's 
standard entropy). This scheme produces a variety of 
ensembles, with q controlling the degree of confinement. 
Some ensembles belong to the Gaussian universality class 
but others exhibit anomalous behavior, characterized by 
distributions having power-law tails. 

The explicit construction of these g-ensembles is pre- 
sented in Sect. ^ where we also derive expressions for 
marginal distributions and the joint density of eigenval- 
ues. Remarkably, for large matrices, the g-ensembles can 
be represented as a superposition of Gaussian ensembles. 
This allows us to obtain closed analytical formulas for 
the eigenvalue density, level-spacing probability distribu- 
tions, etc fSect. lnH) . The comparison of analytical results 
with numerical simulations is the subject Sect. IIVI We 
present in Sect. ^ the concluding remarks. 



II. THE GENERALIZED ENSEMBLES 

For simplicity we restrict our analysis to ensembles of 
real and symmetric matrices H - extensions are straight- 
forward. The volume element in this space is 



N 



N 



(5) 



i=l 



i<3 



where it is understood that matrices are of size N x N. 
Generalized ensembles are obtained by maximizing the 
entropy of Eq. Q subjected to normalization, Eq. 
and 



/ dU trH^ [F(H)]« 
JdU [P(H)]9 



(6) 



with a a constant having units of energy (we are as- 
suming that of H is a Haniiltonian) . Equation © is 
the generalization of the usual constraint that leads to 
the Gaussian ensembles in the standard maximum en- 
tropy approach. Arguments justifying the use of the es- 
cort probabilities P'', and applications of this general- 
ized maximum entropy scheme to various problems can 
be found in Ref. 

Using the Lagrange multiplier technique it is straight- 
forward to find the distribution of maximum entropy: 



P(H) cx exp^ (-A trH^) , 
where we have defined the ^-exponential function pi 



with 



exp (x) = {[1 + (1 - q)x] + } 



&+ — max{&, 0} 



1/(1-9) 



(7) 



(8) 



(9) 



[note that expi{x) = exp(a;)]. The omitted normalization 
constant in Q and the parameter A can be determined 



from the constraints 1^ and (Some preliminary re- 
sults along these lines have been obtained by Evans and 
Michael 

The ensemble defined by Eq. iQ will be called the "q- 
Orthogonal Ensemble" (qOE), as it can be seen in Q 
that the probability distribution depends only on trH^, 
an orthogonal invariant. When q ^ 1 the g-exponential 
function tends to the usual exponential, and one recov- 
ers the Gaussian Orthogonal Ensemble (GOE). Except 
for the q = 1 case, the g-exponential in |(7J) cannot be 
factorized into a product of (marginal) distributions for 
individual matrix elements Hij, which are then corre- 
lated. We can already verify that the cases q < 1 and 
q > I are qualitatively different. Equations © and 
show that for q < 1 the distributions have compact sup- 
port; if (7 > 1, there are always power-law tails (we are 
assuming A > 0, see below). 

To proceed with the analysis of qOE it will be conve- 
nient to think of matrices H as points in a d-dimensional 
euclidean space 

B Ellll- The first N components of 
a point r correspond to diagonal elements Ha, the last 
ones to the upper triangle Hij, i < j: 



(Hi 



,Hnn,V2Hi2,---,V2Hn-i,n) ■ (10) 



The dimension of this space equals the number of inde- 
pendent matrix elements of H, i.e.. 



d 



N{N + 1) 



(11) 



The scaling of Hij by makes the probability distri- 
bution O spherically symmetric in R'^, i.e., PqOE(i") is 
the product of a uniform distribution in the angles, and 
a radial distribution T?| 



P(r;g,a,iV)cxr'*-iexpJ-Ar2) , 



where 



= r • r = trH-^ 



(12) 



(13) 



The observations above imply that qOE belong to the 
wider category of "spherical ensembles" recently studied 
by Le Caer and Delannay 15, 16] . 

For q > 1, the distribution H12|) has a power-law tail 
that goes like l/A^+^^'> with 



A* = 



-d. 



(14) 



Then the normalization condition cannot be satisfied for 
all values of q, but only by those making /i > 0, i.e., 

2 



oo < g < 1 



(15) 



(Note the formal similarity between this problem and the 
generalized random walker in d dimensions [T^.') 
The Lagrange multiplier A is given by 



A = 4 



(t2 2-d{q-l) 



(16) 
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Inside the region (|15|l (i.e., normalizability) A is always 
positive. 

Integrating Eq. lO over all variables but one, we obtain 
the marginal distributions for diagonal and off-diagonal 
matrix elements, 18J 



and 



P{Hii) cx exp^, 
P{Hij) cx expg, 



where 



and 



<? = 



(rf^3)(g-l) 



A' 



d 2-(d-l)(g-l) 



2cr2 



2-d{q-l) 



(17) 
(18) 



(19) 



(20) 



The following properties can be easily verified. The pa- 
rameter q' is an increasing function of and around the 
critical value g = 1 one has 



q + 0[{q-lf] 



(21) 



In addition. A' is always positive. Then, in parallel with 
the global -P(H), the marginal distributions also decay 
as power laws or have compact support, depending on q 
being larger or smaller than one, respectively. We remark 
that the matrix elements are not independent, so, i^(H) 
cannot be reconstructed from the marginal probabilities 

mm- ... 

The joint density of eigenvalues can be obtained in a 
straightforward way: 0, 



N 



N 



P(£i, • • • ,eAr) (X \Y -£»|expJ -A^£- I . 

. . (22) 

The part that is responsible for level repulsion is identi- 
cal to that in GOE because it arises only from orthog- 
onal symmetry. The difference is in the confinement 
term, which in the present case is a non-separable q- 
exponential. Thus the "potential" that confines the spec- 
trum is not a single-particle quadratic well, as in GOE. 
It is rather a mean field, proportional to the moment of 
inertia ^ ef . 

We can get a clear view of the generalized ensembles 
by noting that these are connected to the so-called Fixed 
Trace Ensembles (FTE) and, for large A^, to the Gaussian 
ensembles. In fact, recall that FTE are defined by 0,[ll 



PFTE(H;r, N) cx ,5(trH2 



(23) 



Let /(H) be an arbitrary function and consider the av- 
erages in both ensembles qOE and FTE, namely. 



(/(H))qOE(<Z,a,7V) = y dHPqOE(H;q,a,7V)/(H) , 



(24) 



(/(H))FTE(r,7V) y dHPFTE(H;r,7V)/(H). 



(25) 



Then we have the relation 



(/(H))qOE(<Z,^,A^) = 

dr V(r- q, a, 7V)(/(H))FTE(r, N) . (26) 



The average over qOE can be calculated in two stages. 
First do the average over the angles, for a fixed radius r. 
This correspond to a FTE average. Then average over 
radii, with the weighting function V{r). Of course, the 
same is true for the GOE, which corresponds to the par- 
ticular case q = 1. The relationship between qOE (or 
GOE) and FTE is analogous to that between the canon- 
ical and microcanonical ensembles of Statistical Mechan- 
ics. 

Equation H26|l involves no approximations. Although 
exact, it is not very useful because it requires the knowl- 
edge of fixed-trace averages. However, if one is interested 
in the limit of large matrices, important simplifications 
can be made. 



III. THE LARGE N LIMIT 

The key point is that, for N large enough, the FTE 
average in the r.h.s. of (|26|l can be approximated by 
an average in a GOE having the property (trH^) = r^. 
Then, if we know the GOE average of a given function, its 
corresponding qOE average can in principle be calculated 
by doing just one integration. We will analyze in detail 
two spectral statistics: the eigenvalue density. 



N 



p{e;q,cr,N) ^ (^S{e - 



and the distribution of level spacings. 



p{s;q,a,N) = {S{e. 



i+l 



5)> 



(27) 



(28) 



In the last equation and £,;+i arc two consecutive eigen- 
values lying at the center of the band, i.e., £i ~ 0. It 
is (or will become) clear that other statistics, e.g., two- 
point correlation functions, can also be considered along 
the same lines. 

In order to obtain the qOE averages of H27|) and l|28|) 
we need the corresponding FTE expressions, to be fur- 
ther averaged with 7'(r; q, a, N), as indicated by Eq. H26|l . 
However, we will approximate FTE averages by the corre- 
sponding GOE ones. Then, the basic ingredients become 
the "semicircle law" (for the eigenvalue density), 



p{e; N, r) 



27rr2 



4j,2 



(29) 
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and Wigner's surmise, 
p{s;N,r) = 



O 2 "^^P 



47^7-2 



(30) 



giving the level-spacing distribution. Equations H29|) and 
(|30|l are good approximations for both GOE and FTE 
distributions when N is large [l5l |. 

We recall that if g > 1, the normalization condition 
p5|l limits the value of to a finite domain. On the other 
hand, the case q < 1 does not present such a problem. 
So, we analyze both cases separately. 



A. Ensembles with q < 1 

Except for providing an energy scale, a plays no special 
role. From now on, without loss of generality, we set 
(T = 1. If desired, a can be restored at any time by 
dimensional analysis. 

When N ^ oo {q fixed) the radial distribution of qOE 
tends to 



V{r-q,N) cx r'^-^ [l - r 



„21 1/(1-9) 



(31) 



limited to the domain < r < 1. As c? grows, the distri- 
bution is squeezed against r = 1, being concentrated in 
a smaU region below r = 1, of width 0{l/d) = ©(l/TV^). 
It can be verified that both the level density H29|l and the 
spacing distribution H3U|I . when considered as functions 
of r, have widths which are 0(1). Thus the radial distri- 
bution is much narrower and we can safely approximate 



V{r;q< l,iV->cx)) ~5{r~l) 



(32) 



We conclude that, when q < \ and oo the en- 

sembles qOE tend to the GOE [as far as it concerns the 
distributions being studied, namely Eqs. l(?7)l and (^5)1 ]. 



B. Ensembles with q > 1 

When q> 1 the possible ensembles are restricted to a 
region in the plane q — N that gets thinner as iV — s- oo 
(see Fig. . The natural coordinates in this region are 
N (or d) and [see Eq. (|14|l ]. the latter controlling the 
tails of ^^{r) and other distributions. For instance, sub- 
stituting ((Tl|l and into ((T7|) or ((THl) , one immediately 
verifies that the marginal distributions behave asymptot- 
ically like 



1 



(33) 



The radial distribution H12() . as a function of r, /i, N be- 
comes: 



P(r, /i, TV) cx r 



d-l 



1 + - r' 



(34) 




FIG. 1; When g > 1 the ensembles qOE lie in the region lim- 
ited by the axes and the curve /i = (normalization frontier). 
As A'^ becomes large, the maximum q allowed tends to one. 
Lines correspond to families of ensembles having the same 
power-law tails (labeled by /i). 



This expression allows the identification of some well 
known ensembles as special members of the qOE class: 
the Cauchy-Lorentz ensemble corresponds to fx = 1. An 
integer /i > 1 produces Student's ensembles (see 0, ll^ 
and references therein; see also [20)'l. 

Now we analyze the limit N oo while keeping /i > 
fixed, i.e., we move upwards along the curves of Fig. ^ 
As in the case q < 1, examined before, there is a limiting 
distribution. Some simple algebra leads to 



Vir, fi, N ^oo) cxr~(i+^) exp(-^) 



(35) 



Only when fi ^ oo, V tends to the delta function H32|l . 
and GOE is recovered. For finite /i the width of V^r) 
is at least 0(1). In any case, the average of a given 
GOE distribution with V{r) gives the corresponding qOE 
distribution (via Eq. (|2ni) with (/(H))_ft£; (/(H))GOiJ, 
when A^ oo). Let us first consider the density of states. 
Inserting Eqs. and P5|l into wc obtain 



p[e; m) cx / dr — exp 

J\/Ne/2 ^ 



2r2 



(36) 



This integral can not be expressed in terms of elementary 
functions. However, some information can be extracted 
analytically. Setting e = one obtains the qOE density 
of states at the center of the band. 



jv3/2r[(^ + i)/2] 
r [/x/2] 



(37) 



The behavior for large e can be easily recognized by mak- 
ing the change of variables 2r = \fNez in l|36(l . which 
leads to 



p(e;M) cxe"(i+'') 



dz 



2^1+3 



■ exp 



2^ 



(38) 
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Evidently the tails vanish like e~''^^^\ This is also 
the behavior observed by Cizeau and Bouchaud in their 
"Levy ensembles" of matrices having independent entries 
distributed according to the same law of Eq. ^ T?]. 
We note, however, that the analogies can not be pushed 
further because our ensembles are rotationally invariant 
and Levy ensembles are not (the ensembles of Ref. [Tll | 
belong to the so-called a-symmetric class 0|). 

The calculation of the spacing distribution proceeds as 
before. We have to insert Eqs. ((20)) and into (j^ . 
The result is 



p(s; /i) oc s / drr ''^+'^) 



exp 



2r2 



27r/i 



(39) 

The dependence on s can be easily isolated by a change 
of variables, so we can write: 



p{s;fi) cx s 1 + 



27r/i 



-{l+^L/2) 



or alternatively 



p{s; fi) (X s exp^^ (— as^) 



where 



and 



7V3 



47r /i 



(40) 



(41) 



(42) 



The function of Ea. H4U|l (or Ea. H41|l ') is identical in shape 
with the exact level-spacing distribution of the 2x2 qOE 
having the same /i (see appendix). Then, both distribu- 
tions can be collapsed by a simple scaling of the argu- 
ments. This curious result constitutes a generalization of 
Wigner's surmise to qOE. 

Remark. When analyzing spectral statistics it is usual 
to normalize energies so that the (local) average spacing 
is one (the spectrum is "unfolded"). This amounts to 
measuring energies in units of 



A = 



ds sp{s) . 



(43) 



Note, however, that in qOE the first moment of p(s) does 
not exist for < 1. In these cases, instead of A, one may 
alternatively use the energy scale 



ds s ^ p{s) 



(44) 



Due to level repulsion there is no singularity at s = 0, and 
A always exists, thus representing a characteristic energy 
of qOE. It is close to the inverse of the level density at 
e = 0: 



7rp(0;/i) 



(45) 



with p(0; /i) given in (|37|) . 



IV. NUMERICAL RESULTS 

When thought of as clouds in R'^, via the map of 
Eq. 1)101) . both ensembles qOE and GOE are spherically 
symmetric. This means that qOE can be constructed 
just by rescahng the radii of all points in the GOE cloud 
[iSl iW^ . Thus, the construction of a qOE matrix Hi 
(with parameters q,a,N) can be done in three steps, (i) 
Construct a GOE matrix Ho of size N x N. In this case 
matrix elements are independent and can be calculated 
using Eqs. ((T7I) and ifT^ with q — I. The radius of Hq is 



(46) 



(ii) Choose a radius ri randomly according to the radial 
probability distribution V{ri,q,a,N) of Eq. H12I) . (iii) 
Define Hi as 



Hi 



Ho^ 

To 



(47) 



This is the recipe we followed for constructing qOE matri- 
ces. (If, instead of being a random variable, ri is fixed, we 
obtain a matrix belonging to PTE.) The only difficulty is 
to devise the random number generator, especially when 
V{r) has very long tails. For this purpose we used a com- 
bination of the rejection method and the transformation 
method^ as explained in Ref. |2H . 

In Fig. 13 we show histograms representing densities of 
states obtained from diagonalization of qOE matrices. It 
is clear that they are very well described by the formula 
(|38|l . which was evaluated by direct numerical integra- 
tion. The statistics of level spacings is exhibited in 




FIG. 2: Density of states in the ensembles qOE (normalized to 
one). We compare histograms generated numerically (dots) 
with the theoretical result of Eq. I|88|l (curves). Each his- 
togram was generated from a set of W matrices. We used 
the following values of p: 0.5, 1.5, 2.5, 6.0, oo. Densities with 
larger /x's have larger values at e = and decay faster. In all 
cases A'^ = 20 and a = 1. The dashed line corresponds to the 
GOE semicircle (A'^ oo). 
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Fig. 13 Histograms were obtained by binning data from 
10^ matrices. Each matrix contributed with the "central" 
spacing between levels eiq/2 and 8^/2+1- The analyti- 
cal curves are the g-distributions of Eq. (|40l) . Again the 
agreement between theory and simulations is satisfactory. 
In both figures we observe some small deviations, which 
may be attributed to the relatively small size of the ma- 
trices considered (iV=20). The values of ji were chosen in 
accordance to the following criterion. When ji — 0.5 all 
integer moments diverge. For /i — 1.5 (/i — 2.5) the first 
(second) moment exists but higher ones diverge. The 
case /i — 6.0 is intended to represent an ensemble qOE 
approaching GOE. 



0.8 I ' ' ■ ^ 




s 



FIG. 3: Level-spacing distribution in the ensembles qOE. We 
compare histograms generated numerically (dots) with the 
theoretical result of Eq. (1401 (curves). Each histogram was 
generated from a set of 10^ matrices, each matrix contributing 
one pair of levels. We used the same values of ^ as in Fig.0 
Distributions with larger /i's have higher maxima and decay 
faster. In all cases TV = 20 and a = 1. 



q = 1 giving normal diffusion 0, 0| . 

In the present case, the entropic index q controls 
the confinement, allowing to access different universal- 
ity classes. The anomalies we found can be associated to 
an effect of soft- (or weak-) confinement. However, the 
confining potential in qOE is many-body, as opposed to 
the more common single-particle confinement in standard 
Random Matrix Theories. 

Gaussian ensembles are ergodic |2^. In Sect. lIIII we ex- 
plicitly used the fact that individual large qOE matrices 
have GOE statistics. This implies that qOE is not ergodic 
(if fi is finite). In random matrix theory, non-ergodicity 
is considered to be a drawback because, it is argued, pre- 
dictions (ensemble-averages) are compared with data ob- 
tained from a single system. We do not object to this 
reasoning, but just mention that in some cases empirical 
data are indeed extracted from ensembles of Hamiltoni- 

ans mm. 
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APPENDIX A: GENERALIZED WIGNER 
SURMISE 

Here we apply to the level spacing distribution in 
qOE with N ^ 2. We want to calculate 

p{s; q, a) = {S [(£2 - £1) (H) - s])qOE , (Al) 



V. CONCLUSIONS 

This paper explores the possibility of using a modi- 
fied maximum entropy method to construct random ma- 
trix ensembles. Whereas the usual logarithmic entropy 
leads to ensembles characterized by exponential laws, the 
power-law entropy of Eq. naturally produces ensem- 
bles with long-tailed distributions, (e.g., the so-called 
g-Gaussians [13 )■ These ensembles (gOE) recover the 
Cauchy-Lorentz and Student's ones for special choices of 
the parameters. 

Of course, the same families of ensembles can also be 
obtained from say the standard maximum entropy ap- 
proach, but at the expense of introducing a complicated 
constraint. A similar situation arises in the maximum 
entropy approach to anomalous diffusion, where one can 
choose between an unappealing constraint and a 
nonextensive entropy [iTl Il9j | . By sweeping the param- 
eter q one switches from sub- to superdiffusive regimes. 



where ei < £2 are the eigenvalues of H. First we need 
p{s;r) for FTE. This is a known result [T6| : 



p(s;r) 



1 



V2r V2r2 - 



(A2) 



(if the argument of the root is positive; zero otherwise). 
According to (|^ . p(s; q, a) is obtained by averaging (|A2|I 
with the radial weight of qOE, Eq. l(T^ with d = 3. We 
show the result for the case q > 1: 



where 



and 



p{s;q,a) cx sexp^,, (-A"s^) , 



1 + 9 



A" 



3-9 ' 



3-g 



40-2(5/3 _q) 



(A3) 



(A4) 



(A5) 
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Equation (|A3|) plays the role of Wigner's surmise for In Sect. HTll this expression is compared with the spacing 
qOE. It is more useful to rewrite HA3() in terms of the distribution for the large N case (Eq. (|40ll '). 
parameter /x of Eq. l(T^ . The result is 

/ 3s2\-(1+m/2) 

p(s;g,a)(xs 1 + — . (A6) 
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